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ABSTRACT 



Context. A previous analysis of methanol and formaldehyde towards the Orion Bar concluded that the two molecular species may 
trace different physical components, methanol the clumpy material, and formaldehyde the interclump medium. 
Aims. To verify this hypothesis, we performed multi-line mapping observations of the two molecules to study their spatial distribu- 
tions. 

Methods. The observations were performed with the IRAM-30m telescope at 218 and 241 GHz, with an angular resolution of ~ W" . 
Additional data for H2CO from the Plateau de Bure array are also discussed. The data were analysed using an LVG approach. 
Results. Both molecules are detected in our single-dish data. Our data show that CH3OH peaks towards the clumps of the Bar, but 
its intensity decreases below the detection threshold in the interclump material. When averaging over a large region of the interclump 
medium, the strongest CH3OH line is detected with a peak intensity of ~ 0.06 K. Formaldehyde also peaks on the clumps, but it is 
also detected in the interclump gas. 

Conclusions. We verified that the weak intensity of CH3OH in the interclump medium is not caused by the diff'erent excitation con- 
ditions of the interclump material, but reflects a decrease in the column density of methanol. The abundance of CH3OH relative to 
H2CO decreases by at least one order of magnitude from the dense clumps to the interclump medium. 

Key words. ISM: individual objects (Orion Bar) - ISM: abundances - ISM: molecules - ISM: structure 



1. Introduction 

Given its proximity and nearly edge-on orientation, the so- 
called Orion Bar is one of the clearest examples of a photon- 
dominated region (FDR). For this reason, the Orion Bar has 
been extensively observed in past decades to test theoretical 
models of FDR structure, chemistry, and energetics. Its nearly 

■ edge-on orientation make direct observations of the gas strat- 
ification possible, from ionised gas to neutral atomic gas to 
molecular gas as a func tion of the increasing distance from the 
ionisation source (e.g.. iTielens et al.l Il993l:1van der Werf et al.l 

; 119961: Ivan der Wiel et al.l l2QQ9 l). These studies show that the 

■ molecular gas consists of clumpy molecular cores embedded 
in an interclum p gas. While the clumps have de nsities of sev- 
eral 10^ cm-3 (lYoung Owl et al.1 12QQ0I: |Lis_&Schilke 2003), 

■ the interclump material has lower densities (10"^- 10^ cm"^, 
' iHogerheiide et al.lll995l:lYoung Owl et"aD l2000'). Such observa- 
tional results_^e_succes^ by theoretical works 
(e.g. jGorti & Hollenbachll2002h . 

Observations of methanol (CH3OH) and formaldehyde 
(H2CO) in the Orion Bar have been r eported in the pas t 
('Hogerheiide et al."l995l: lJansen et al.1 1 9951: iLeurini et al.l2 006). 
^Leurini et al. ( 2006) studied the excitation of both molecular 
species in the line of sight of one molecular clump through 
a multi-line analysis at 290 GHz. Their findings (high density 



and relatively low temperature for CH3OH, high temperature, 
and relatively low density for H2CO) suggest that methanol and 
formaldehyde do not trace the same material, but the first is 
found in the dense gas associated with the clumps, the second in 
the warmer and less dense gas of the interclump material. Since 
both molecular spe cies can efl&ciently form on grain surfaces 
(iHidaka et al.l 12004 ). the authors suggested that photodissocia- 
tion of methanol to form formaldehyde (Le Teuflf et al. 2000) 
takes place in the interclump medium, while the high density 
shields CH3OH in the clumps and prevents its photodissocia- 
tion. Whereas grain surface reactions are the only viable route to 
methanol, formaldehyde can also form in the gas phase via the 
reaction CHs-fO (e.g., Le Teuflf et al. 2000), which could con- 
tribute to the interclump abundance of H2CO. 



Send offprint requests to: S. Leurini 

* Based on observations carried out with the IRAM-30m tele- 
scope and the Plateau de Bure interferometer. IRAM is supported by 
INSU/CNRS (France), MPG (Germany) and IGN (Spain). 



However, the study of 'Leuri ni et al.l (l2006l) on CH3OH and 
H2CO is based on low angular resolution (~ 20''), single- 
pointing data. To confirm that the two species indeed trace diflTer- 
ent media, we mapped a large area of the Orion Bar in formalde- 
hyde and methanol at 218 and 241 GHz, respectively, with an 
angular resolution of ~ IV\ We also analysed existing interfer- 
ometric data for the 21 8 GHz formaldehyde lines, to compare the 
observed fluxes from single-dish and interferometer and derive 
constraints on the size of the emitting region. This work rep- 
resents a first observational attempt to study the distribution of 
CH3OH and H2CO in the Orion Bar and to verify their diflTerent 
origin. 
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Fig. 1. In grey scale, the integrated intensity of the H2CO (3o,3 - 2o,2) line (left) and of the CH3OH-A (5o - 4o) transition (right) 
observed with the IRAM-30m telescope. Black contours (4, 8, and 12 tim es 0.032 Jy beam"^ km s"^) show the distribution of the 
H^^CN (1-0) emission observed with the Plateau de Bure interferometer (iLis & Schilkdl20Q3l) . For H2CO the contours start from 
1.8 K km s"^ in steps of 1.8 (equal to 3(t), for CH3OH from 0.3 K km s"^ in steps of 0.3 (equal to 3cr). In the left panel, the dashed 
lines indicate the strips used for Fig.O the solid lines outline the regions used to derive the average H2CO and CH3OH emission 
in the interclump material (Figs. |3(b)| and |4(c)| ); the dotted lines mark the region shown in the right panel. In the right panel, the 
numbers label the positions of the clumps detected in H^^CN; the dotted circles outline the mosaic of seven fields observed with the 
PdBI (seeiO. 



2. Observations 

2.1. IRAM-30m telescope 

The observations were performed in Fe bruary and March 20 07 
using the HERA multi-beam receiver dSchuster et al.l l2004l) at 
the IRAM-30m telescope. The HERAl and HERA2 pixels were 
tuned to 241.850 GHz in lower side band (LSB), to detect the 
CH3OH (5;^ - 4k) band. On February 23 and 24, the HERA2 re- 
ceivers were tuned to 218.349 GHz (LSB) to detect the H2CO 
OKa,K, - '^Ka,Kc-i) band. The backend used was VESPA with 
a bandwidth of 160 MHz and a resolution of 0.3125 MHz for 
the CH3OH setup (corresponding to a velocity resolution of ~ 
0.4 km s"^), 160 MHz and 1.25 MHz for the H2CO setup (corre- 
sponding to ~ 1.7 km s"^). In this way, the CH3OH (5k-4k) band 
is fully covered with the only exception being the ^ = - 1 , 0-E 
lines. The H2CO setup covers the 3o,3 - 2o,2 and 32,2 - 22,1 lines, 
while the 32, 1 -22,o transition lies outside the observed frequency 
range. 

We used the on-the-fly mode with a rotation of the multi- 
beam system of 9.7°, to ensure a Nyquist sampling between 
the rows. The reference position used as the centre of the map 
was Qf20oo = 05^35^^25^3,^2000 = -05°24'34'^0, correspond- 
ing to the "Orion Bar (HCN)" position of Schilke et al. (2001), 
the most massive clump seen in H^^CN (Lis & Schilke 2003), as 
well as the target of the spectral survey of .Leurini et al. ( 2006). 
The OFF position was chosen to be (500'', 0'') from the centre 
of the map. 

The pointing was checked on several continuum sources and 
found to be accurate within 2''. The observations were done 
under varying weather conditions, and the Tsys span the range 
540-1200 K for HERAl and 700-1600 K for HERA2. The half- 
power beam width of the IRAM-30m telescope at 218 GHz is 
~ IT' and ~ 10'' at 242 GHz. We used a beam efficiency of 



0.55 at 218 GHz, and of 0.50 at 242 GHz to convert antenna 
temperatures T* into main-beam temperatures TmtEl- 

The region covered by our observations is shown in Fig.[T] 

2.2. Plateau de Bure interferometer 

We observed the H2CO 3o,3 - 2o,2 line at 218.2 GHz with the 
Plateau de Bure interferometer in March, April, and December 
2004. We observed a 7-field mosaic centred on CK2000 = 
05^35^^25.280^ ^2000 = -05°24'39".300. The field positions 
followed a compact hexagonal pattern to ensure Nyquist sam- 
pling in all directions. The imaged field-of-view is almost a cir- 
cle with a radius of 22.5". On March 28 and 29, the six antennas 
were used in the 6 Cp configuration with 3 and 1.5 mm precip- 
itable water vapour, which translated into system temperatures 
of 450 K and 250 K. On April 22, the six antennas were used 
in the more compact 6Dp configuration with 6 to 10 mm PWV, 
leading to a system temperature of about 1000 K. On December 
12, the six antennas were used in the 6Cp configuration, with 
3.5 mm PWV, leading to a system temperature of about 700 K. 
Taking the time for calibration and data filtering into account, 
this translates into an on-source integration time of useful data 
of 5.3 hours for a full 6-antenna array. The typical 1 mm resolu- 
tion is 2.4". We used the 30m data described above to produce 
the missing short spacings. 

The data processing was done with the GILDAS0 software 
suite (Pety 2005). Standard calibration methods implemented 
in the GILDAS/CLIC program were applied using close bright 
quasars (0528-hl34 and 0607-157) as phase and amplitude cal- 
ibrators. The absolute flux was calibrated using simultaneous 



^ http://www.iram.es/IRAMES/telescope.html 
^ See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILDAS softwares. 
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Fig. 2. Top: uv coverage for one field of the hybrid 30m+PdBI 
mosaic. Bottom: Associated dirty beam. Please note that the 
dirty beam has many secondary side-lobes as high as 40% of 
the main lobe. 



measurements of the PdBI primary flux calibrator MWC349. uv 
tables were then produced at a relatively coarse velocity resolu- 
tion of 1.7 km/s. 

All other processing was performed with the 
GILDAS/MAPPING software. The single-dish map from 
the IRAM-30m was used to create the short-spacing pseudo- 
visibilities unsampled b y the Pla teau de Bure interferome- 
ter (Rodriguez-Fernande z et al.ll2Q08) . These were then merged 
with the interferometric observations. Each mosaic field was 
imaged independently, and a dirty mosaic was built thr ough a 
linear combination of these dirty images (l Guethetal.i ri995). 
The dirty mosaic was then deconvolved using the standard 
Hogbom CLEAN algorithm. While all these procedures are 
standard and usually easy to use in the GILDAS/MAPPING 
software, we had to take special precautions with this data 
set because the short usable integration time at PdBI implied 

1) a low signal-to-noise ratio of the interferometric data and 

2) dirty beams for each field with large secondary side-lobes 
(see Fig. [2l). As the 30m data show that the extended signal 
fills most of the field-of-view observed with the interferom- 
eter, we subtracted the spectrum averaged over the observed 
field-of-view from the 30m data before any processing and we 
added this averaged spectrum again after the deconvolution 
of the 30m-hPdBI hybrid data set to recover the correct flux 
scale. This simplifies the deconvolution considerably by the 
usual CLEAN algorithms (i.e. much fewer CLEAN components 
are needed) because it avoids the deconvolution of extended 
uniform intensities. Second, we tapered the visibility weights 



Table 1. Detected methanol and formaldehyde transitions in the 
IRAM-30m telescope data 



line 


y 

[MhlzJ 




clump 1 


interclump 


H2CO 3o,3 - 2o,2 


218222 


21.0 


ya 


Y 


H2CO 32,2 ~ 2i 1 


218475 


68.0 


ya 


Y 


CH3OH-A 5o- 4o 


241791 


34.8 


ya 


Y 


CH30H-A5±3 -4±3 


241832 


84.6 


yb 


N 


CH3OH-A 52-42 


241842^ 


72.5 


yt 


N 


CH30H-£ 53-43 


241843 


82.5 


yb 


N 


CH30H-£5i -4i 


241879 


55.9 


ya 


N 


CHsOH-E 5+2-4+2 


241904 


~ 60^ 


yb 


N 



^ detected at the original resolution of the IRAM-30m telescope 
data 

^ detected at a resolution of ~ 20" 

blended with the 53 - 43-E 
^ Eup =60.7 K for the 5_2 - 4_2 line, 57.1 K for the 5+2 - 4+2 tran- 
sition. 



with an axis- symmetrical Gaussian of 80m FWHM to get a more 
symmetrical beam (going from A.V'xl.T' for natural weighting 
to 3.3'' X L8'' after tapering) and to improve the signal-to-noise 
ratio for the extended structures. The final noise rms measured 
at the centred of the mosaic is about 0.3 K in channels of L7 
km/s width. This leads to a maximum signal-to-noise ratio of 18 
for the hybrid data cube. However, the produced data cube has 
its brightness dynamical range limited by the poor dirty beam. 
The cleaned data cube was finally scaled from Jy/beam to 
temperature scale using the synthesised beam size. 

3. Observational results 

3.1. single-dish data 

Figures [3(a)l and [4(a)l show the spectra at 218.4 and 241.8 GHz 
towards the central position of the map. Both formaldehyde tran- 
sitions are detected. Only the 5o - 4o-A and 5i - Ai-E methanol 
transitions are detected at the original resolution of ~ 10''. 
However, when increasing the signal-to-noise ratio by smooth- 
ing to a lower resolution of ~ 20", CH3OH transitions wi th up - 
per level energies lower than 85 K are also detected (Fig. |4(b)| ). 
The detected lines of both molecular species are reported in 
Tableffl 

The region mapped with the IRAM-30m telescope is pre- 
sented in Fig. [TJ where in the left panel we show the integrated 
intensity of the strongest detected H2CO line, while the right 
panel shows the emission from the strongest CH3OH transition 
along the Orion Bar alone. For comparison, the integ rated inten- 
sity of the H^^CN (1-0) from iLis & Schilkd (12003 1 is overlaid 
on our data. Both molecules clearly trace the molecular medium 
associated with the Orion Bar, and peak to the southwest of it. 
This peak probably belongs to the low-intensity bridge that con- 
nects Orion South to the Orion Bar. This structure was detected 
in previous observations at velocities between 7 and 9 km s"^ 
(e.g., Tauber et al. 1995). The southwest peak of CH3OH and 
H2CO detected in our data has a peak velocity of 8 km s"^. In 
addition, formaldehyde is detected towards the north of the map, 
at a posit ion that spatial ly coincides with continuum emission at 
350 yum (iLis et al .11 19981) and with a velocity between 9 and 1 1 
km s"^ 

Along the Bar, both molecules trace the elongated region 
where molecular clumps are located, with an absolute emis- 
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Fig. 3. Spectrum of the H2CO Ok,,k, 
from the IRAM-30m telescope. 
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Rest Frequency (MHz) 

(b) 

'^Ka,Kc-i) band towards clump 1 (top) and averaged over the interclump medium (bottom) 



sion peak on clump 3 of ' Lis & Schilkd (l2003h . Moreover, both 
species are detected on a secondary peak, northeast of clump 1 
(0^2000 = Q5^35^29. 6^(52onn = - ( )5°23[ 59^^9), not covered by 
the observations of Lis & Schilkd (12003 V but close to a peak of 
CO(6-5) (Fig. O, which could represent an additional clump of 
dense gas along the Orion Bar. To better visualise the overall 
morphology of the Or ion Bar, in Fig. [5] we show the 20 cm con- 
tinuum emission from Yusef-Zadel which traces the ion- 
isation front; the CO(6-5) integrated emission (Lis e t al .11 1998!) . 
which shows the temperature distribution of the molecular gas; 
and the H^^CN (1-0) emission, which traces the dense clumps. 

In Fig.[3]and|4lwe compare the H2CO and CH3OH spectrum 
towards clump 1 to the corresponding spectra averaged over the 
region of the interclump gas outlined in Fig. [T] (left). The inte- 
grated intensity of the CH3OH (5o - 4o)-A transition in the inter- 
clump medium is a factor ~ 9 lower than in the clump, the one 
of H2CO 32,2 - 22,1 line a factor of 7 (see Table [2l). For a better 
visualisation of the formaldehyde and methanol morphologies, 
we show in Fig. [6l the variation in the integrated intensities of 
the CH3OH-A (5o - 4o) and H2CO (32,2 - 22,1) lines for three 
strips in the map; the exact location of th e strips is shown in 
Fig.[T] The strip across clump 1 (Fig. |6(a)| ) suggests a larger ex- 
tension of the H2CO emission with respect to CH3OH, while 



the variation in the integrated intensity of the two species along 
the Bar is very similar (Fig. |6(c)| ). Although we present data for 
the weakest of the two H2CO lines, the signal-to-noise ratio in 
the methanol data is still lower than for H2CO and could bias 
our results. From these strips and from the comparison between 
the H2CO emission and the CO(6-5) map, it also emerges that 
formaldehyde extends to a greater distance from the clumps in 
the molecular cloud than towards the ionisation front. We also 
smoothed the data to lower resolutions, to verify whether the 
CH3OH and H2CO emission are similarly affected by beam di- 
lution. At a resolution of 50'', the H2CO emission is still ~ 37% 
of the intensity of the original data, while at the same resolution 
the CH3OH peak intensity drops down to ~ 18% of the original 
value. These tests suggest that formaldehyde emission is associ- 
ated with the clumps but also with the interclump medium, while 
methanol emission is only found in the dense molecular clumps. 
In Sect, m we investigate whether the different morphologies of 
the two molecular species stem from an excitational or observa- 
tional bias, or whether they imply different abundances of H2CO 
and CH3OH in the two media. 

Finally, the spectral resolution used for the methanol ob- 
servations allows us to study the velocity field along the Bar. 
The channel maps of the (5o - 4o)-A line (Fig. [71 smoothed 
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Fig. 4. Spectrum of the CH3OH (5y^ - 4^^) band towards clump 1 at a resolution of ~ 10 (top), smoothed to a resolution of - 
(middle), and averaged over the interclump medium (bottom). Data are from the IRAM-30m telescope. 



20' 



to 0.77 km s ^) show that clump 1 peaks arou nd 9.5 k m s ^ 
and clump 3 around 10.5, as for H^^CN (Lis & Schilki[2003). 
The secondary peak detected to the northeast of clump 1 has a 
peak velocity around 11 km s"\ while the molecular cloud to 
the south of the Bar is red-shifted with respect to the clumps 
(visr ~ 8 km s"^). This trend in velocity (blue-shifted velocities 



in the northeast, red- shifted values in the southwest) was also 
found by Young Owl et al.. (.2000) in their HCN maps. 

3.2. Hybrid interferometric+single-dish data 

Figure [8] compares the 30m (top panels) and hybrid 30m-hPdBI 
(bottom panels) data cubes. The left column displays the images 
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Table 2. Line parameters 
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J TuBdv 












[K] 


[km s-i] 


[km s-i] 


[K] 


[Kkms-i] 


[km s-i] 






H2CO (3o,3 - 


2o,2) 




CH3OH (5o - 


4o-A) 


clump 1 


3.2 


8.9 


2.0 


1.1 


1.8 


1.4 


interclump 


0.3 


1.2 


3.1 


0.06 


0.2 


4.0 



^ corrected for the spectral resolution, Av = ^(Avobs)^ - (Avres)^ 




R.A. 



(a) 

Fig. 5. Distribution of the CO(6-5) peak brightness temperature 
(colour image). The black contours (from 40% of the peak inten- 
sity, in steps of 10%) show the 20 cm continuum emission, which 
traces the ionisation front. The dotted contours (from 15% of the 
peak intensity, in steps of 10%) are the integrated intensity of the 
H2CO 3o,3 - 2o,2 line; the blue contours (as in Fig. [T]) represent 
the dense clumps; the square marks the position of the secondary 
peak identified in H2CO and CH3OH. 

integrated between 8.3 and 11.7 km/s. The middle column dis- 
plays the spectra averaged over the whole interferometric field- 
of-view. The hybrid spectra are identical to the 30m ones within 
the noise level, confirming that all the extended emission filtered 
out by the PdBI is correctly recovered from the single-dish data. 
The right column of Fig.[8]displays the spectra averaged over the 
central clump marked in the left column. The hybrid spectrum is 
brighter than the 30m spectrum implying beam dilution in the 
single-dish data. 

From the PdBI data and the hybrid 30m-l-PdBI data, we can 
therefore conclude that the H2CO emission comes from both an 
extended component and a compact one, still unresolved at the 
resolution of the IRAM-30m telescope. 

4. Discussion 

4.1. Formaldehyde 

Line ratios of formaldehyde lines can be used to infer the physics 
of th e gas in dense molecular regions ( Mangum & Wootten 
Il993h . Because of the coarse spectral resolution of our observa- 
tions (Av ~ 1.7 km s"^), the line profiles of the H2CO transitions 
are poorly resolved. As a result, we do not have a correct mea- 
sure of the peak intensity of the lines, diluted over the beam, but 
our values are a lower limit to the true main-beam line intensi- 
ties. However, the areas are conserved quantities and can be used 
to study the excitation of H2CO in the Orion Bar. 



Figure [9] shows the ratio of the integrated intensities of the 
3o,3 - 2o,2 fine to the 32,2 - 22,1 fine. The analysis is limited to 
the region where the signal-to-noise ratio of the 32,2 - 22,1 inte- 
grated intensity is at least equal to 3. The ratio increases from 
the northeast, where the clumps are, to the southwest. 

We ran LVG calculations for column densities of H2CO be- 
tween 10^^ and 10^^ cm"^, temperatures in the range 20-200 K, 
and densities of lO'^-lO^ cm"^, and calculated the behaviour of 
the;^^ value of the observed 3o,3 - 2o,2 and 32,2 - 22,1 fine inten- 
sities as function of these quantities. Two possible regimes are 
found (see Fig. [T0|, independent of temperature. At low densi- 
ties (n< 10^ cm"^), the;^^ value depends only on the product of 
the H2 density and H2CO column density (low-density branch). 
For higher densities, the value varies only slowly with H2 
density and is mostly a function of H2CO column density (high- 
density branch). 

For the values in the interclump gas, assuming a beam fill- 
ing factor of unity, we find a nominal minimum for the in- 
terclump medium at T = 180 K, n(H2) = 6.6 x 10^ cm'^ 
A^H2C0 = 1.4 X 10^^ cm"^. However, as can be seen in Fig. [TOl 
these parameters are not well constrained. Using the Icr confi- 
dence level, we derived a lower limit to the kinetic temperature 
of the interclump equal to 76 K, and the product x A^H2C0 is 
equal to ~ 3 x 10^^ cm"^ for dens ities smaller than 10^ cm"^. By 
comparison with PDR models, i Young Owl et al.l (I2QO O) found 
densities of h ydrogen nuclei of lO '^-lO^ cm"^ in the interclump; 
similarly, Hog erheijde et al.l (Il995h found an interclump H2 den- 
sity of 3 X 10"^ cm"^ by modelling mm line observations. Thus, 
by using a molecular hydrogen density between 5 x 10^ and 
5x10^ cm"^ we can constrain the column density of H2CO to 
values in the range 6x10^^-6x10^"^ cm"^. The high-density 
branch does not appear to be relevant for the interclump medium. 

In Sect. I3.2l we analysed the combined PdBI-h30m data and 
concluded that the H2CO emission is composed of two compo- 
nents, one extended and one unresolved at the resolution of the 
30m telescope at this frequency. However from these data, it is 
difficult to infer the exact size of clump 1 and the contribution 
of the interclump gas to the total emission of H2CO. In Fig. O 
we marked clump 1 with an ellipse of 15'' x 5'' (P.A.=47°) and 
showed that the IRAM-30m data are aff'ected by beam dilution at 
this scale, implying a smaller size for the clump. Therefore, for 
clump 1 , we corrected the main-beam line intensities for the con- 
tribution from the interclump and for the beam dilution assum - 
ing a typical size for the clump of 7'' (see iLis & Schilkell2003h . 
As expected given the similar value of the line ratio on clump 1 
and in the interclump, the values derived are similar to those of 
the interclump: T = 170 K, n(H2) = 6.6 x 10^ cm-^ A^H2C0 = 
1.6 X 10^"^ cm"^. The Icr contour gives a value of 3.8 x 10^^ cm"^ 
for the product x A^H2C0 for A^H2C0 > 5 x 10^^ cm"^ and 
^(H2) < 5 X 10^ cm"^. For the low-density branch, we derive 
a lower limit of 50 K for the kinetic temperature. For the high- 
density branch, more appropriate for the clump, the column den- 
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Fig. 8. Top line: IRAM-30m data. Bottom line: hybrid 30m+PdBI deconvolved data cube. Left column: Images of the integrated line 
intensity of the H2CO (3o,3 - 2o,2) transition between 8.3 and 11.7 km/s. The values of the contour levels are shown on the colour 
scale. Middle column: Spectra averaged over the whole interferometric field-of-view marked as the red circle on the images. Right 
column: Spectra averaged over the central clump marked with the blue ellipse on the images. 



sity is constrained to values 2.8 x 10^"^ - 5.4 x 10^"^ cm"^. A lower 
limit of 80 K can be derived for the kinetic temperature. It may 
appear surprising that the temperature is so little constrained, 
although the rat io of these lin es has been used as a tempera- 
ture tracer (Man g um & Woott en 1993; Hogerheijde et al. 1995). 
However, for the values of the line ratio we observe, the LVG 
predictions show a more complex behaviour with temperature 
(see also Fig. 13 of lMangum & WoottenlfT993h . Taking the mea- 
surement errors into account by calculating the value, as we 
do, and using the line intensities and not just the line ratios then 
constrains only a lower limit to the temperature. 

The results of the LVG analysis of H2CO are presented in 
Tables 

4.2. Methanol 

iLeurini et aP (l2QQ4l) have studied the excitation of the CH3OH 
{5k - 4^) band as function of the temperature, density, and 
methanol column density of the gas for a range of physical con- 
ditions applicable to the Orion Bar. They concluded that line ra- 
tios within this band are not sensitive to the temperature of the 
gas for temperatures higher than 30 K, but only depend on the 
density. 

The methanol emission from cl umps 1 and 3 is studied in 
a separate paper (iParise et al. | [2Q09>) : clump 1 is found to be 
warmer than clump 3, although the errors are large (3(T ranges: 
Ti ~ 45^^!^ K, 73 ~ 35^1^ K), but with similar column densities 
and densities (A^cHbOH ~ 3 x 10^"^ cm"^, ^ > 5 x 10^ cm"^). These 
values were derived by modelling the 5y^ - Ak and 6k- 5k bands 
observed with the IRAM-30m (observations presented in this pa- 
per) and the APEX telescopes, respectively. For the models, the 
IRAM data were smoothed to the resolution of the APEX data 
{FWHM htdim ~ 2C). A source size of IC was used, since sev- 



eral clumps fall in the beam of the observations. This is equiv- 
alent to assuming that all clumps have similar physical condi- 
tions. Since the CH3OH emission was found to be optically thin, 
we can correct the results for a source size of 1'' (as assumed for 
the model of H2CO) and derive the equivalent column density. 
This corresponds to a total column density of 5.4 x 10^"^ cm"^, 
or to 2.7 X 10^^ cm"^ for CH3OH-A, assuming that the two sym- 
metric states have the same column density. The Icr confidence 
level is 2.3 x lO^^ - 3.1 x lO^^ cm'^. 

For the interclump medium, we analysed the CH3OH spec- 
trum averaged over the area of the interclump medium out- 
lined in Fig. O We ran models for the excitation of methanol 
in the range of densities 10^ - 10^ cm"^, column densities 
IO12 _ iQis and temperatures 20 - 200 K. The observed 

main-beam brightness temperature of the (5o - 4o)-A line in the 
interclump medium is 0.06 + 0.01 K. Assuming that the emission 
comes from an extended source (e.g., beam filling factor ric = 1), 
we find a situation similar to the one described for H2CO and we 
can identify two diff'erent regimes. In the first (n < 10^ cm"^), a 
fit to the data is found for Ncn^on-A < 1.3 x 10^"^ cm"^: the col- 
umn density decreases with increasing density, but their product 
is almost constant (between 2 x 10^^ cm"^ 50 K < T < 200 K, 
and 2 x 10^^ cm'^ 20 K < T < 50 K). For higher densities, 
the 5o - 4o line thermalises. In this regime, the column density 
increases slowly (A^cHgOH-A = 1.3 x 10^^ - 2 x 10^^ cm"^). For 
both regimes, the dependence on the kinetic temperature is not 
strong and the results presented are valid for temperatures in the 
range 20-200 K. However, increasing the temperature of the gas 
implies a decrease in the column density of methanol for the 
low-density case (T > 50 K, Nch^qu-a < 3 x 10^^ cm"^, see 
Fig.in]). 

As already discussed in the previous section, from the litera- 
ture we know that the density of the interclump gas is less than 
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Table 3. Summary of the model results 





A^H2C0-p 




[CH30H-A]/[H2CO-/?] 






[10^3 cm-2] 




clump 1 


(28 - 54)^ 


23-31 


0.4-1.1 


interclump 


(6 - 60)^ 


(0.04 - 4)^ 


7 X 10-4 - 0.7 






(0.04 - 0.4)^ 


(7 X 10-4 - 7 X IQ-^y 



^ assuming n = 10^ cm~^ 
^ assuming « = (5 X 10^ - 5 x 10^) cm-^ 
assuming T > 50 K 



10^ cm ^, and therefore that the low-density regime applies to 
the modelling of the interclump. The best fit to the data is found 
for A^cH30H-A ~ 3 X 10^2 cm'^, T ~ 80 K, ^ ~ 9 x 10^ cm'^ For 
densities in the range 5x10^ - 5x10^ cm"^, and for the values 
of the product x ^ given above, we derive column densities in 
the range 4 x 10^^ - 4 x 10^-^ cm"^. However, if we allow only 
temperatures above 50 K, as found in our data and by other au- 
thors, the inferred methanol column density in the interclump is 
4 X 10^1 < < 4 X 10^2 cm-2 (see Table©. 



5. Abundances of methanol and formaldehyde in 
the Orion Bar 

From the analysis performed in the previous section, we can 
compute the abundance of methanol relative to formaldehyde in 
the interclump medium and in the dense clumps, and verify that 
the drop of intensity of the methanol emission in the interclump 
medium is the result of the diff'erent physics of the interclump 
relative to the dense clumps, or whether it reflects a real decrease 
in the CH3OH abundance. 

Given an estimate of H2 column density, we could also infer 
the abundance of both molecules relative to molecular hydrogen. 
An estimate of the H2 column density can be derived from the 
continuum emission at (sub)millimetre wavelengths, assuming a 
given temperature for the dust and optically thin emission. The 
continuuni emis sion at 350 yum of the Orion Bar was studied by 
iLis et al.l (1 19981) . These authors defined an average temperature 
for the dust in the OMC-1 regi on of 55 K on the basis of its 
FIR colours (IWerner et al.| ll976') and excluded colder tempera- 
tures for the Orion Bar given the good agreement of the 350 yum 
continuum emission and the CO(6-5 ) emission, wh i ch ori ginated 
in the outer PDR layers. However, I Werner et al.l (11976^ found 
an average FIR colour temperature of 75 K in the Orion Bar, 
which could imply that the Orion Bar is significantly warmer 
than the rest of OMC-1. For a dust temperature of 55 K in the 
interclump medium, and for the mean value of the 350 fim emis- 
sion over the region used to derive the average spectra of the 
interclump medium (Fig. [T]), we derive a column density of H2 
of 5.3 X 10^^ cm"^. This translates into to 3.4 x 10^^ cm"^ for 
Td = 75 K. A value of Q- 101 g cm"^ was used for the dust opacity 
(lOssenkopf & Henningll 19941) . For clump 1, for a temperature of 
45 K (as derived from methanol), and under the assumption that 
the dust and the gas are thermally coupled at the high density of 
the clumps (e.g., Krti gel & Walmsleyll 19841) . the H2 column den- 
sity is 2 X 10^^ cm"^. This is not the average column density on 
the beam, but it was corrected for a source size of 7'' and for the 
contribution of the interclump medium at the position of clump 
1. These values translate into abundances of H2CO in the range 
(11 - 110) X 10"^ for the interclump,(1.4 - 2.7) x 10"^ for the 
clumps. For methanol, we obtain an abundance relative to H2 of 



(0.08 - 8) X 10"^ for the interclump, (1.2 - 1.6) x 10"^ for the 
clumps. 

However, these estimates are aff'ected by large uncertainties: 
a diff'erence of 20 K in the dust temperature implies an uncer- 
tainty of a factor of ~ 2 in the estimate of the H2 column den- 
sity in the interclump, uncertainty that could be even greater if 
the average temperature in the region of the interclump used 
in our analysis is higher than 75 K. Similarly, the estimate of 
the H2 column density on the clumps is aff'ected by large un- 
certainties because of the assumption that the dust is thermally 
coupled to the gas at high densities and because of the large er- 
rors on the temperature of the gas from the methanol spectrum 
(28 < Ti < 92 K, 20 < < 52 K). 

On the other hand, the uncertainties in the determination 
of the abundance of CH3OH relative to H2CO are related only 
to the errors in the estimate of the column density of the two 
molecules, which already include the uncertainties in the tem- 
perature. Therefore, we believe [CH3OH - A]/[H2C0 - p] to 
be more solidly estimated than [CH3OH - A]/[H2] or [H2CO - 
/^]/[H2]. Using the column densities derived in the previous sec- 
tion, we infer an abundance of CH3OH relative to H2CO of the 
order 0.4-1.1 in the dense clumps, and 7x10""^- 7x10"^ in the 
interclump. If we assume that the temperature of the interclump 
medium is higher than 50 K, then abundance of methanol drops 
down to 7x10""^- 7x10"^. Since it is a reasonable assumption 
that the temperature of the interclump medium is 50 K or even 
higher, we conclude that the abundance of methanol relative to 
formaldehyde decreases by at least one order of magnitude in 
the interclump medium in comparison to the dense clumps. 

A possible explanation for the decrease in abundance of 
methanol relative to formaldehyde from the clumps to the inter- 
clump medium is photodissociation in a low-density enviroment. 
Although the photodissociation ra tes of CH3OH and H 2CO are 
of the same order of magnitude (iLe Teuff et al.ll200Ql) . one of 
the products of the photodissociation of methanol is formalde- 
hyde. Therefore, both molecules can be destroyed by the ra- 
diation field in the interclump medium, but H2CO may be re- 
plenished through the photodissociation of CH3OH. However, 
large uncertainties aff'ect the measurement of the photodisso- 
ciation rates, and so more accurate models would be required 
to test this hypothesis. Alternatively, as already discussed by 
Leurini e t al.l (|2006), H2CO could be formed in the interclump 
medium through gas phas e reactions, whic h are not efficient for 
the formation of CH3OH (iLuca et al.ll2QQ2b . 

6. Conclusion 

In a previous analysis of methanol and formaldehyde towards 
the Orion Bar, we suggested that the two molecules might trace 
diff'erent environments and that, while CH3OH is associated with 
the clumpy molecular cores of the Bar, H2CO is found in the in- 
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Fig. 7. Channel maps of the CH3OH-A (5o - 4o) emission be- 
tween 11.3 and 7.4 km s"^ observed with the IRAM-30m tele- 
scope. 
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Fig. 6. Variation in the integrated intensity of the CH3OH-A (5o- ^ig. 9. Ratio of the integrated intensity of the H2CO 3o 3 - 2o 2 
4o) and H2CO 32,2 - 22,1 Imes observed with the IRAM-30m transition to the 322 - 221. The solid black contours show 



telescope across clumps 1 (a), the interclump medium between 
clump 1 and 2 (b) and along the Bar (c). In Fig. [T] we outline the 
strips used to extract the plots. 



the integrated intensity of the 3o,3 - 2o,2 line (levels are from 
1.8 K km s"^ in steps of 1.8). Typical errors along the Bar are of 
the order of 0.3. 



terclump material. To test this hypothesis, we mapped the Orion 
Bar in both molecular species with the IRAM-30m telescope, 
with a spatial resolution slightly higher than the expected size of 
the clumps. Additional data were taken with the IRAM Plateau 
de Bure Interferometer in H2CO. 

Both molecules are detected in the Orion Bar in our single- 
dish data. Our data show that CH3OH peaks towards the clumps 
of the Bar, but its intensity decreases below the detection thresh- 



old in the interclump at individual positions. By averaging over 
a large region of the interclump medium, the strongest of the 
CH3OH lines in our setup (5o - 4o-A) is detected with a peak 
intensity of ~ 0.06 K. On the other hand, contrary to the hypoth- 
esis formulated in our previous study, formaldehyde is detected 
towards the clumps and the interclump gas. 

Using an LVG program, we studied the excitation of H2CO 
and CH3OH in the Orion Bar. We suggest that formaldehyde is 
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Fig. 10. square distribution of the 3o,3 - 2o,2 and 82,2 - 22,1 
line intensities in the [^(H2), A/^(H2C0)] plane for a temperature 
of 100 K for the two sets of observations clump and interclump 
shown in Fig.O The plot shows, for each spatial region, the low- 
density branch and the high-density branch discuss ed in the text. 
Note that to derive the column densities given in § 14.11 one has 
to multiply with the line widths. 
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Fig. 11. Results of the statistical equilibrium calculations for 
CH3OH. The solid and dashed lines show the line intensity of 
the 5o - 4o-A line as measured towards the interclump medium 
(0.06 K) for temperatures of 20 and 80 K, respectively, as func- 
tions of density and column density of CH3OH-A . 



present in both components of the Bar (clumps and interclump 
material), with column densities up to 5.4 x 10^"^ cm"^ in the 
clumps and between 6x10^^-6x10^^ cm"^ in the interclump. 
These values only refer to para-formaldehyde. From the analy- 
sis of methanol, we concluded that the reason for the drop of 
intensity of CH3OH in the interclump medium is not the dif- 
ferent physical conditions with respect to the clumps, but a real 
drop in its column density compared to the clumps (down to 
4 X 10^^ cm"^ for temperatures above 50 K). The abundance of 
methanol relative to formaldehyde decreases by at least one or- 
der of magnitude in the interclump medium compared to that in 
the dense clumps. 



Despite the large errors on our estimates, our observations 
reveal that the column density of methanol and formaldehyde in 
the clumps are of the same order of magnitude, while the col- 
umn density of methanol in the interclump is lower than that 
of formaldehyde. This may be a result of photodissociation of 
CH3OH in the unshielded interclump gas, or it may reflect that 
H2CO can be produced in the gas phase more efficiently than 
CH3OH. Definitive observational conclusions would require a 
much deeper integration towards the interclump gas and a bet- 
ter characterisation of its density and temperature. More detailed 
chemical models of PDRs, including grain surface reactions, are 
clearly needed to properly interpret these observational results. 
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